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We report adjustable magnetic 'bouncing' and focusing of a dilute 87 Rb Bose gas. Both the con- 



densate production and manipulation are realised using a particularly straight-forward apparatus. 
The bouncing region is comprised of approximately concentric ellipsoidal magnetic equipotentials 
with a centre that can be adjusted vertically. We extend, and discuss the limitations of, simple 
Thomas-Fermi and Monte-Carlo theoretical models for the bouncing, which at present find close 
agreement with the condensate's evolution. Very strong focusing has been inferred and the obser- 
vation of atomic matter-wave diffraction should be possible. Prospects look bright for applications 
, in matter-wave atom-optics, due to the very smooth nature of the mirror. 

o . 

^ ■ I. INTRODUCTION 

The realisation of gaseous atomic Bose-Einstein condensation (BEC) in 1995 1] sparked intense theoretical and 
experimental interest in this hitherto inaccessible quantum degenerate regime. Several landmark experiments have 
followed, and activity in the field continues to escalate Q. There has been a similar heightened interest in the area 
of experimental atom-optics, following the demonstration of a variety of techniques for manipulating cold atoms 

As BECs are the ultimate coherent cold atom source, it has become important to develop high-quality atom- 
optical elements for manoeuvring them without destroying their coherence. Potentials typically used for atom-optical 
manipulation arise from magnetic and/or optical dipole forces. We will focus in particular on the use of cold atoms in 
conjunction with atomic mirrors and lenses. Flat £j and curved @ optical dipole mirrors, as well as flat 0,0,111 an d 
curved [|| magnetic mirrors have all been realised with conventional cold atoms in static confi gura tions. D yna mic 
optical dipole mirrors lQj, flat 0,0] and curved dynamic magnetic mirrors, as well as ID and 3D jl5l IH| 
' dynamic magnetic lensing have all been demonstrated. The temporal variability provided by such mirrors and lenses 
\Q . adds an extra dimension of control. However, dynamic magnetic mirrors and lenses tend to be softer than their static 
counterparts - i.e. atoms have a longer interaction distance with the mirror potential. The optical quality of the 
original static mirrors has been refined and a natural extension is the use of these relatively hard mirrors as 

atom-optical elements to focus and image coherent atomic sources such as BECs. 

First steps in this direction were made in the experiment of Bongs et al. |17| , where the ingredients of a BEC and 
atomic mirror were successfully combined with the demonstration of condensate reflection from a flat spatial dipole 
mirror based on a rapidly rastered laser beam. Our mirror has the additional advantages that it can be inelastic, 
with an adaptable radius of curvature. Here we give a more detailed description of our earlier work |l8l Il9| on the 
reflection and focusing of a BEC from this soft curved temporal magnetic mirror. A similar mirror has been used 
to reflect, focus and store atom laser beams 20]. The ability to adaptively focus and launch condensates will enable 
the realisation of BECs with very high densities, and it has already enabled us to create extremely cold and coherent 
\ collimated atomic clouds. 

An important aspect of our mirror design is that all magnetic elements are a relatively large distance (d > 20 mm) 
from the atoms. Thus, although the magnetic fields are relatively weak compared to hard atomic mirrors in which the 
atoms make a close (d < 100 /jm approach to the magnetic surface, microscopic corrugations in the mirror potential 
are drastically reduced because curvature in the magnetic field scales as d~ 3 . Such large-scale mirrors thus have the 
advantage that they are extremely smooth, ensuring fully coherent atomic manipulation. The roughness of our mirror 
was measured to be less than 1 mrad rms |l9j , where this upper bound was probably only limited by the sensitivity 
of the measurement technique, and our atomic mirror is therefore one of the smoothest in the world [P9ll2l| . 

We begin the discussion with an overview of our 'simple' BEC apparatus, followed by benchmarking experiments 
that were performed on the condensates. We will then concentrate on bouncing and focusing BECs: contrasting 
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and discussing the limitations of two theories, paying particular attention to the theory of experimentally attainable 
situations. 



II. 'SIMPLE' BEC CREATION AND BENCHMARKING 

In the Sussex BEC experiment N > 10 9 87 Rb atoms can be collected at the low pressure end of a low-power double 
magneto-optical trap (MOT) |22j system (Fig. QJ. The three extended-cavity diode lasers providing the necessary 
780 nm (5s 2 Si/2 — > 5p 2 P 3 / 2 ) MOT light have a very simple, stable design and can be cheaply manufactured from 
commercially available components with minimal machining |23| . A single diode laser originally provided all trap 
laser light for the BEC experiment (11 mW total per MOT in 25mm diameter (FWHM=20mm) beams), however 
two independent lasers are now used for the high pressure and low pressure MOTs - increasing the low pressure MOT 
atom number with a minor cost to simplicity. For both MOTs the trap laser light is detuned A = — 13 MHz (i.e. red) 
of the F = 2 — * F' = 3 transition, whilst the repumping light is resonant with the F = 1 — > F' = 2 line. 

Every 80 ms the « 10 7 atoms that have accumulated in the high pressure MOT are transferred to the low pressure 
MOT simply using the light pressure from a 150 /is push pulse of saturated light resonant with the F — 2 — > F' = 3 
transition. The shape of the apparatus (Fig.^) and the relatively small angular spread of the transferred atoms allow 
the realisation of a moderate (rj rs 40%) transfer efficiency without requiring a magnetic guiding mechanism. The low 
pressure MOT is loaded with Nlp ~ 10 9 atoms in 30 s, however faster loading times could be achieved by inserting 
an aperture in the transfer tube - increasing the present 5 : 1 MOT pressure ratio. 




FIG. 1: The experimental setup, denoting the high (HP) and low (LP) pressure MOT chambers. Vertical laser beams are not 
shown. The 20L/s ion pump is all that is needed to maintain the vacuum. 

We find that a dark/compressed MOT phase [24| is unnecessary for BEC production, however this small complication 
would no doubt improve the condensate population. The atoms in the low pressure MOT are cooled to 40 /xK in optical 
molasses (A = —35 MHz), then optically pumped into the \F, mp) = 12, 2) state using a circularly-polarised light pulse 
resonant with the F = 2 — » F' = 2 transition. A Ioffe-Pritchard (IP) |2J| magnetic trap is subsequently formed around 
the atoms. The constant magnetic field in the LP MOT chamber can be altered in < 1 ms using three orthogonal 
pairs of magnetic coils - thus providing the necessary field for molasses and optical pumping as well as positioning 
the LP MOT centre relative to the IP trap. 

Absorption imaging of the atoms is performed with the optical pumping laser beam. The imaging lens is an / = 5 cm 
doublet designed to minimise distortion, with a magnification of M — 0.80(1). The short (8 /is) pulse of near-resonant 
imaging light propagates horizontally at an angle of 9 w ir/6 to the — x axis of the magnetic trap. We refer to the 
resulting horizontal axis in our images as —x' (Fig.^l. The imaging/optical pumping beam, like the push beam, is 
derived from the high pressure MOT trap laser. The intensities of all laser beams derived from the two trap diode 
lasers can be altered in times t < 1 /xs via their corresponding constant-frequency acousto-optic modulators (AOMs). 
Complete extinction of all laser light (important during magnetic trapping) is realised using mechanical shutters. As 



3 



the free-running rms laser linewidth is Az/ = 300 kHz (for times t < 200 ms), the laser frequencies can be precisely 
adjusted over ±300 MHz simply by unlocking the laser and making small, abrupt changes to the laser current. 

The Ioffe-Pritchard magnetic trap coils are formed from water-cooled 3 mm o.d. copper tubing: a 9-turn 'cuboidal' 
baseball coil with average side dimensions 45 x 45 x 55 mm 3 , and two three-turn 64 x 101 mm 2 'rectangular' bias coils. 
The coils' 10 V, 250A constant-voltage power supply was controlled by two simple but effective integrator circuits 
driving independent water-cooled banks of MOSFETs [2^|. Series rms current noise levels of AI/I < 5 x 10~ 5 were 
obtained, and the absence of shunt resistors across both the baseball and bias coils ensured the maximum level of 
frequency- insensitive magnetic noise cancellation when the coils were operated completely in series |27f. The magnetic 
trap's residual bias field of 10 G could be lowered to 0— 10 G using a pair of 'extra' bias coils. With a typical compressed 
magnetic bias field of 1 G the axial and radial trapping frequencies were measured to be v z = 10 Hz and v r = 223 Hz 
respectively at a magnetic trap current of 220 A. 

BEC experiments commonly use two or more large-throughput vacuum pumps, however due to the relatively small 
surface area of our apparatus the vacuum is essentially maintained by a single 20L/s ion pump (Fig. The lifetime 
of the magnetic trap, r = 70s, is certainly adequate for obtaining condensates. 

After a 32 s RF evaporative cooling ramp pure condensates are reproducibly obtained with more than 10 5 atoms. 
Two 'bench-marking' experiments were performed on our BEC: the temperature-dependent grou nd-state occupation 
was compared with the results of Refs. [2j| EM |3(], E| and the BEC's ballistic expansion |32L l33| was contrasted 
with a Thomas-Fermi [34J theoretical model j3f|. Our data yielded good agreement with these earlier results, and 
more details of our ground-state occupation experiment will be given shortly. For our ballistic expansion experiments, 
absorption imaging was only valid for condensate ex pan sion times t > 10 ms (cf. |23) due to the expected lensing 
effects 36]. As in Ref. 33] . non-destructive imaging 36] will be used in future to probe the condensate's behaviour 
at short release times. 

A variety of models have been used by different groups to model the condensate and non-condensate fractions of 
cold atomic clouds. Three density distributions commonly used in BEC data analysis will be summarised here: 
• the Gaussian distribution 
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The density distribution obtained from absorption imaging is integrated along the line of sight of the probe laser 
beam. By using the co-ordinate transformation into the x'y'z' (imaging) frame 
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and integrating along the absorption imaging beam axis (a;'), one obtains the equivalent 'two-dimensional' density for 
each of these distributions: 
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where the width parameters oy, a z , , and the maximum TF radii A y , , A z , have the same angular dependence, i.e. 

cr^sin 2 ^, A y , = A y , A z , = J A z 2 cos 2 8 + A x 2 sin 2 8 . (5) 



(jy, = <7 y , a z , = \l a z 2 cos 2 ' 
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In the Sussex BEC experiment the viewing angle of the camera is 9 w ir/6 because the probe laser optical access to 
the vacuum cell is constrained by the IP magnetic trap coils. In order to obtain a viewing angle along one of the 
magnetic trap axes it would be necessary to directly overlay the probe laser beam with the MOT laser beams. This 
could be done in future by placing polarising beamsplitters in front of the MOT quarter-wave plates - thus allowing 
light to be coupled into the MOT beam path, then coupled out again after passing through the vacuum cell. 

The different models used by previous groups to measure the condensate ground state occupation used either a 
bi-modal Gaussian distribution 30, 31], a 'Bose'/Thomas-Fermi distribution for the thermal/condensate populations 
|29| . or a bimodal Gaussian distribution with a separate Gaussian fit to the 'wings' of the distribution to measure 
the temperature psj . We found that, for our trap parameters and within experimental uncertainties, the atomic 
distributions were fit equally well with a Gaussian/Thomas-Fermi or Bose/Thomas- Fermi distribution. We found 
that bi-modal Gaussian fits were consistently worse compared to these two models. This was expected due to the 
high condensate density (no(t = 0) « 2 x 10 14 cm 3 ) in our magnetic trap, which ensured the applicability of the 
Thomas-Fermi approximation. Fig. illustrates the Gaussian/Thomas-Fermi and Bose/Thomas- Fermi fits. In the 
left graph the temperature is derived directly from the Gaussian or Bose distribution fitted to the thermal atoms, 
whereas in the right graph the temperature is determined by fitting a Gaussian to the 'wings' of the distribution (like 
in Ref. 0). 
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FIG. 2: Relative condensate population as a function of temperature. The data was obtained by fitting bi-modal distributions 
to the average of live imaged thermal/condensate clouds after a 23 ms ballistic expansion. Data points of increasing temperature 
correspond to final RF evaporative cooling frequencies of 710, 715, 760 kHz. Thomas-Fermi distributions were fit to the BEC 
clouds. In the left image both the temperature and number of atoms in the thermal cloud was obtained from a fitted Gaussian 
(triangles) or Bose (diamonds) distribution. In the right image the atomic populations were obtained in the same way, however 
the temperature was obtained by fitting a Gaussian to the wings of the distribution. 



III. BOUNCING AND FOCUSING CONDENSATES 



Atoms moving adiabatically in a magnetic field experience a potential, U(r) = — // B = gprnpfisB, that depends 
on the magnitude of the magnetic field, B, as well as the strength and sign of their atomic state's mag netic moment 
in the direction of B (gpmpfiB)- Initial magnetic focusing experiments with neutral atoms |l4Ul5tlla ] have involved 
cylindrically symmetric systems comprised of one or two circular coils. Judicious application of pulses of magnetic 
potential with positive axial curvature (and negative radial curvature) followed by a magnetic potential pulse with 
neg ative axial curvature allows one to three-dimensionally focus a cold cloud of atoms which passes through the coils 

Ideally, one would like to focus the atoms with a single magnetic pulse as the dynamics are much simpler. The 
BEC experiment already relies on a magnetic potential with positive curvature in all three spatial dimensions - i.e. 
the magnetic trap - this trap is therefore ideal for 3D atom focusing, and for large enough fields it can also be used 
as a mirror. Another advantage of using the same coils for BEC formation and magnetic bouncing/focusing is that 
one does not need to precisely position the respective centres of two different coil systems. A simpler, separate mirror 
consisting of two coils in a quadrupole configuration can be built [ 1 .31] , however along the axis of the coils the magnetic 
field magnitude varies linearly, resulting in only 2D focusing. 

A final pair of magnetic coils was used to generate a vertical bias field B c = {0, B c , 0} whilst the atoms are 
'bouncing' on the IP potential. The time-dependent magnetic field to second-order is thus: 



B = So{0,0,l} + Bi{aj,i/ c -»,0} 

+B 2 /2{-xz, -yz, z 2 - (x 2 + y 2 )/2} 



(6) 
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where the IP trap's magnetic bias, gradient and curvature (Bq, B\, Bi respectively) are proportional to the series 
current in the trap coils, and the centre of the magnetic trap is shifted vertically upward by the variable amount 
y c = B c /B\. The magnitude of the magnetic field is parabolic in the axial (z) direction, but hyperbolic in the 
co-ordinate r = x 1 + (y — y c ) 2 and thus linear for radii r > \/2Bq/Bi. 

In the experiment, a condensate is released at time t = and expands ballistically until t = 15 ms when both the 
IP trap (at w 70% of the 220 A maximum current) and a constant vertical bias field (B c = — 400 G) are pulsed 
on for 5 ms. After 30 ms an identical magnetic field pulse allows the realisation of a second 'bounce'. During the 
'bounce' pulses the atoms experience a magnetic field magnitude that is weakly confining (parabolic) in the z direction, 
strongly confining (with strength dependent on B c ) in the x direction and almost linear in the y direction (see Fig. 1 of 
Ref. |l9j)- Linear spatial dependence of the magnetic field magnitude results simply in a change of the centre-of-mass 
motion (bouncing), whereas parabolic dependence yields focusing/defocusing of the atomic cloud (Fig. 

During reflection the COM motion changes slightly as the vertical magnetic bias is altered because its field has a 
small amount of curvature. Horizontal motion arises partly due to a tilt in the IP trap with respect to gravity, and 
also due to the fact that the magnetic potentials have an increasingly large slope in the z direction as y decreases. 
This yz term in the magnetic field magnitude also causes the cloud to begin rotation in the x = plane during the 
bounce. 

Atoms bouncing in a gravity-cavity are unstable when they are dropped from a height greater than half the radius 
of curvature of the bouncing surface (h > R r j2) [3^ |. 




FIG. 3: One dimensional bouncing on a pulsed potential. The left images ((a), (c), (e)) show the trajectories of atoms with zero 
initial velocity and a range of heights. The right images ((b), (d), (f)) show the trajectories of atoms with zero initial position 
and a range of velocities. Images (a) and (b) show reflection from a linear potential. Images (c) and (d) show reflection from a 
parabolic potential if the drop height h is equal to the radius of curvature of the mirror R. Images (e) and (f) show reflection 
from a parabolic potential if the drop height is h — R/2. 



We have experimentally studied situations where the drop height compared to the mirror radius of curvature in 
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the x direction is h = R x and h — R x /5.4 respectively, with the results illustrated in Fig. 2 of Ref. [ijj. Due to 
our imaging angle, the instability of the condensate in the x direction is not seen in the z' zcos(-7r/6) + xsin(7r/6) 
direction until after the second bounce. 

We now turn to the theoretical analysis of the bouncing experiments. It has been shown [35| that a BEC, in the 
Thomas-Fermi regime and experiencing a time-dependent parabolic potential U(r,t) = ^ X^=i (i) 2- * 
atomic spatial distribution which obeys the simple scaling law: 



n(r, 0) = max 



n 



where X ri (t) is the solution of 



d 2 X n (t) 2 w n (0) 



(7) 



dt 2 



+ u> Ti (t) 2 x ri (t) = ° 5 ( 8 ) 



•V; (0) = 1 and the initial BEC maximal radii are A Ti . We note that if the parabolic potential has a time- varying 
centre r c (t), i.e. if U(r,t) = ^ Yli=i u n ( r c(t)) 2 (fi — r ic{t)) 2 , then the same scaling law applies for the BEC if one 
subtracts out the centre-of-mass (COM) motion r (t) and uses the magnetic curvatures acting at the COM. 

The second theoretical model is a Monte-Carlo simulation. Initial positions for j = 1, 2, M atoms distributed on 
a spherical Thomas- Fermi distribution, n(r) oc max[l — r 2 , 0], can be obtained by inverting the equations arising from 
the cumulative density function — 

Pe<jjr oc J d6 J snv<f)d(f) J r 2 n(r) dr — pg p^ p r — 9 cos (r 3 /3 — r 5 /5) 

— i.e. by uniformly choosing 0j G [0, 27r], p rj S [0, 2/15], and p t j }j S [—1,1] and by inverting p r . and p t j }j to obtain rj 
and <pj = arccos(p0 i ) . 

The initial cartesian positions of the atoms rj = ri, vm in a cylindrically symmetric trap (A x = A y = A r ) are 
then given by 



;(A r sin(</)j) cos(8j), A r s'm((j>j) sm(9j), A z cos(</> :) )), 



r 3 = r J 

where A r = ^^/^, A z = -^y^ and p(N) = {2-Kf/ b m 1 / b {lhh 2 Nav r 2 v z ) 2 ^/2 is the chemical potential. The 
BEC's spatial distribution has a td = 15 ms expansion in which the co-ordinates rj are scaled by Eq. |S| with cor- 
responding velocities Vj = {^^(t)xj, ^-(t)yj, ^-{t)zj}. For times t > 15ms, when the density of the condensate 
and thus the effect of interatomic repulsion has dropped dramatically, the atoms are allowed to follow their classical 
ballistic trajectories. 

Our magnetic coils were simplified to line segments, and the Biot-Savart law was then used to generate the magnetic 
field and its Taylor expansion. Wc could thus ensure that in our case no important third- or higher-order terms are 
neglected in the Monte-Carlo model if one uses the magnetic field magnitude from approximation Eq. For the 
Thomas- Fermi model we use Eq. [^Jto derive the local magnetic field curvatures - which vary slightly with the centre- 
of-mass position of the BEC during the bounce pulses. 

In the Monte-Carlo simulation the condensate maintains an approximately Thomas-Fermi distribution (Eq. in 
the yz' imaging plane throughout the bouncing sequence. However, due to the yz term in the magnetic potential, the 
distribution is rotated in the yz plane by an angle <f> yz (or (j) yz i in the yz 1 plane). In order to find the radii, and tilt 
angle (j) yz ', of the Monte-Carlo distribution, the cartesian co-ordinates of the M Monte-Carlo 'atoms' are projected 
into the yz' plane with co-ordinates (j/i, z^). It can be shown that given the moments of the cloud in the yz 1 plane 

M M M -, M , M 

2 



i=l i=l 

one obtains the rotation angle of the cloud: 

<h vz i — — arctan 



(yz') 



?' 2 _ v 2 
■'rms ilrms 



where z' rrns = sj (z' 2 ) — (z') 2 and y rms — \J (y 2 ) — (y) 2 are the rms 'radii' of the cloud in the non-rotated yz' frame. 
Additionally, in the yz' frame rotated by (f> yz > , the true rms radii of the condensate are given by 



V z' rms 2 cos 2 (j) yz > + y r ms 2 sin 2 <f> yz , + (yz') sm(2(/) yz >) 
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Y rms = y z' rms 2 sin 2 (t> yz > + y rms 2 cos 2 tf> yz i - (yz') sin(20 yz /) 

and one can show that these rms radii can be converted to Thomas-Fermi distribution maximum radii by multiplying 
by a factor of y7- 

The theoretical models work well |19| in describing the absorption images taken from our current viewing angle, but 
they also show that even more interesting behaviour should be seen in the xy plane. When B c = OG (i.e. h = R x ) 
then after the first bounce the condensate comes to a tight focus in the x direction. The Thomas- Fermi model (Fig.0J, 
with its inclusion of atomic interactions but neglect of anharmonic spatial magnetic fields, predicts a very tight focus 
with an x cloud maximum radius of A x = 18 nm. The de Broglie wavelength of the BEC at the focus, 38 nm, can be 
obtained from its velocity. Thus, if a suitably harmonic mirror is used and the focused matter wave has a transverse 
width smaller than the matter wavelength, one may expect to see atomic diffraction during the focus . 



8C 



400 






i 

1 

i 


300 




,'\ 
/ \ 


I 


200 




/ \ 
/ 3r" 
/ \ 


\ 

i 


100 




^'V \ 

- -/ \ 

/ \ 


i 

i 

I 



20 40 60 80 

t (ms) 



FIG. 4: The time-evolution of the BEC's x, y and z radii (A x dash-dot-dotted, A y dash-dotted and A z dashed) in microns 
obtained using the Thomas-Fermi model. Also of interest is the scaled relative density hit) = 100 no(i)/no(0) oc l/(A x A y A z ), 
which, following a slow decay, has two sharp peaks when A x briefly comes to a tight focus. The initial peak density in the 
condensate is no(0) — 2 x 10 14 cm -3 , and at a time t{ = 24.2 ms (shortly after the condensate bounces for the first time), 
the peak density is four times larger for a duration of m 10 fis. Subsequent foci of A x are even tighter, however they depend 
strongly on the harmonicity of the mirror. 

We now compare the Thomas Fermi model with three versions of the Monte-Carlo model. The Monte-Carlo models 
use: the full magnetic potential originating from Eq. (MC), the local harmonic approximation to the magnetic 
potential at the centre of mass (MCh), and the harmonic MCh model with the inclusion of the effects of the x 2 y local 
potential aberration term (MCu x 2 y)- We note that the MCh model and Thomas-Fermi model are identical during the 
initial condensate expansion (t < 15 ms), and it is only during focusing - if the density is large enough for interatomic 
repulsion to bc important - that the models may differ. 

In Fig.[S]we look at the behaviour of the foci in the xy plane as predicted by the various Monte Carlo models. We find 
that the MC model has a maximum density of no = 9(2) x 10 14 cm 3 in a 'bow-tie' shaped focus which varies location 
with time. The MCh model has an extremely tight elliptical focus, leading to a peak density of no — 24(2) x 10 14 cm 3 
- this density is more than ten times higher than the Thomas-Fermi model in which interatomic repulsion limits the 
density at the condensate focus. The MC^ x 2 y model has a similar (but crisper) focal shape and a similar maximum 
density (no = 10(2) x 10 14 cm 3 ) to the MC model this illustrates that the 'bow-tie' focus is mainly generated by the 
x 2 y terms in the magnetic potential from Eq. The approximate (FWHM) temporal duration of the above Monte 
Carlo density maxima during the focus is 30, 2, 30 /zs for the MC, MCh, and MG^ x 2 y models respectively. This can 
be compared to the 10 us duration of the focus seen in the Thomas- Fermi model. In addition, the foci of the Monte 
Carlo models occur approximately 10 us before the Thomas-Fermi focus, as interatomic repulsion delays the focal 
time. 

The Thomas- Fermi model thus predicts very tight foci in the xy plane, however the anharmonicities of our present 
mirror mean that our Monte-Carlo model will be better at predicting the shape of the xy focus. The Monte-Carlo 
models' neglection of the interatomic repulsion which is included in the Thomas-Fermi model, also appears to be 
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FIG. 5: A theoretical comparison of the temporal foci of the atomic distribution in the xy plane subsequent to the first bounce. 
Three different Monte-Carlo models are used (as described in the text): MC, MCh and MCn^s^ models are illustrated in the 
upper, middle and lower rows of the figure, respectively. Each image has x x y (horizontal x vertical) widths of 4 /im x 200 (jm, 
and the images correspond to times (left to right) t — tt — 30 /is, t{ — 20 fis, t{ + 10 fis where tf = 24.2 ms is the time at which 
the focus occurs in the Thomas-Fermi model. 

invalid, as the MCh model differs dramatically from the Thomas-Fermi model at the focus. During the focus very 
high densities are obtained, but only for a very short duration compared to the time over which interatomic repulsion 
is important during the initial condensate expansion (w 2 ms as can be seen in Fig. [2J. Thus although we might 
expect the MC models to break down during the focus, the expansion trajectory of the condensate after the focus 
may well be valid. 

When we begin probing such extremely small length scales in the xy plane (nanometres - rather than the > 1 fim 
diffraction-limited resolutions of present BEC imaging systems), both of the simple theories described here - the 
Thomas-Fermi and Monte-Carlo models - will therefore cease to be valid. These experiments will need to be compared 
to calculations using the more complicated time-dependent 3D Gross-Pitaevskii equation, however it is worth noting 
that the Thomas-Fermi model should be adequate if a sufficiently harmonic magnetic mirror is employed, and the 
Monte Carlo model may well suffice for describing condensate behaviour prior to and immediately after focusing if 
the tightness of the focus does not induce atomic matter-wave diffraction. 

In conclusion, we have created a relatively simple and robust BEC apparatus. With this apparatus we have 
demonstrated adjustable focusing and reflection of a BEC from a soft, extremely smooth, and adaptive temporal 
magnetic mirror. If highly parabolic magnetic mirrors of this type are used, this technology should enable one to 
focus atoms down to nanoscales. In order to investigate such focusing, we plan to image the condensates in the xy 
plane in future, with the hope of observing atomic matter-wave diffraction. Highly parabolic mirrors can also create 
extremely cold, collimated atomic beams for use in atomic clocks and we have already demonstrated a rudimentary 
matter- wave 'beam expander' in this vein [l9|. 
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